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A model of superconductivity is proposed taking into account repulsive particle interaction, meso- 
scopic phase separation and softening of crystalline lattice. These features are typical of many high- 
temperature superconductors. The main results obtained for the model are: (i) phase separation 
is possible only if repulsive forces play a significant role; (ii) the critical temperature as a function 
of the superconducting phase fraction can have non-monotonic behaviour; (iii) superconductivity is 
possible in heterophase systems even when it would be forbidden in pure samples. These results are 
in agreement with experiments. 

I. INTRODUCTION 

In the standard theory of superconductivity one employs a reduced Hamiltonian involving only the attractive part 
of an effective particle interaction, responsible for the existence of superconductivity, and omitting the repulsive part 
of the interaction as irrelevant. This approach was adapted in the original paper of Bardeen, Cooper and Schrieffer 
[1]. The consideration of the Coulomb interaction leads to a renormalization of the coupling constant [1,2], which can 
be effectively included in the formula for the superconducting critical temperature [3]. In usual metals the averaged 
Coulomb interaction is small as compared to the coupling parameter mediated by phonon exchange. A detailed 
discussion of reasons why the repulsive interaction may be neglected has been given in a review by Legget [4] . 

These arguments concern the direct influence of repulsive interactions. However, they can affect superconductivity 
indirectly. An obvious example is that all interactions are included in the definition of the electronic energy bands and 
the corresponding density of states [5,6]. An interaction that is not directly responsible for the pairing mechanism, 
can, nevertheless, drastically strengthen the order by providing a high density of states at the Fermi surface, thus 
ensuring a large value of the effective coupling parameter, which, in turn, leads to a higher superconducting transition 
temperature. Such an increase of the density of states can occur, for instance, when a van Hove singularity is located 
almost precisely at the Fermi level [7]. 

It is necessary to keep in mind that it is Coulomb interactions, and not pairing interactions that define the space 
structure of matter. The electrical properties change drastically from one structure to another [5-8]. Moreover, the 
stability or instability of a crystalline (and chemical) structure is closely connected with superconductivity. This was 
noticed long before the discovery of the high - temperature superconductors. For example, Testardi [9] analysed a 
large number of superconducting A-15 compounds and concluded that superconductivity and structural instability are 
really related: the higher the superconducting critical temperature, the more probable is the occurrence of structural 
instability. Indeed, most superconductors with Tc ^ 20K clearly exhibit instability. 

The observation that higher superconducting Tc is almost always related to increased instability received strong 
confirmation after the historic discovery of high - Tc superconductors by Bednorz and Miiller [10] . Possibly, the most 
careful, critical, material systematics for high - Tc superconductors has been given by Phillips [11]. He stresses that if 
we review the known high - Tc superconducting materials we are unavoidably confronted with mechanical instabilities 
as a factor which always accompanies Tc , whether in intermetallic or oxide - based superconductors. This occurs so 
often that we would hard put to dismiss the repeated appearance of such instabilities as mere coincidence. According 
to Phillips [11], there is no room for doubt that it is lattice instabilities (and nothing else) that produce high - Tc 
superconductivity just as much in cuprates as in intermetallic compounds. 

The main result of structural instabilities is an essential softening of the lattice. This softening is manifested 
in various precursor effects occurring near the superconducting transition temperature, such as anomalies in elastic 
moduli, in strain parameters, and in sound velocity [9,11]. The softening of phonon modes can be observed by 
infrared, Raman, and neutron scattering experiments [11]. Nuclear gamma - resonance studies [12,13] also show that 
a dramatic lattice softening occurs at Tc for high - temperature superconductors, the most noticable result of which 
is an anomalous sagging of the Mossbauer factor. 

The review of experimental data also provides evidence that the lattice softening and the increase in Tc are both 
associated with the formation of phase precursor clusters with structural (spatial or chemical) composition different 
from the basic component [9,11]; the linear size of the clusters is of the order 10 — lOOA . Such mixed structures 
exist in many ceramic superconductors, as has been observed by using pulsed-neutron scattering, synchrotron x - 
ray powder diffraction, nuclear magnetic resonance, nuclear quadrupole resonance, and nuclear gamma resonance 
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[5,11,14]. For instance, the low - temperature orthorhombic and high-temperature tetragonal structures coexist in 
oxide superconductors in a region around [15-17]. A more detailed discussion of these and other experiments can 
be found in a recent review [18]. The relation between the existence of soft vibrational modes and the appearance of 
clusters of a mesoscopic size has been also proved by molecular dynamics computer simulation for metastable glassy 
models [19,20]. 

The clusters of one phase structure inside another, occurring near a phase transition point, are randomly distributed 
in space, and often, being unstable, change with time. To emphasize the randomness of their distribution, they can be 
called hetero structural, or heterophase fluctuations, and to stress their mesoscopic sizes, they can be named mesoscopic 

structural fluctuations [18]. 

Since structural characteristics are intimately related to electronic and conducting properties, the existence of 
structural fluctuations implies the appearance of spatial fluctuations in the superconducting properties. Mesoscopic 
structural fluctuations are accompanied by large energy fluctuations which favour a stochastic separation of a sample 
into superconducting and normal regions [18,21]. There is clear experimental evidence that the superconducting 
cuprates display a nanoscale phase separation into insulating and superconducting nanodomains [7]. This is in 
agreement with some simple models [7,22-24] predicting that such a phase separation can be thermodynamically 
profitable. Recent experiments with high - temperature superconductors confirm that only part of a given sample - 
often less than a few percent - is in a superconducting phase (see [25-27] and references therein). 

In the present paper we suggest a way of taking into account the three mutually interrelated factors: Coulomb 
interaction, phase separation, and lattice softening; all of which factors are very important in high - temperature 
superconductors. To clearly emphasize the influence of these factors as such we shall use the standard approximations 
accepted for superconductors. We give a detailed analysis of the superconducting critical temperature as a function 
of parameters related to the attracting and repulsive interactions and to the superconducting phase fraction. 

Everywhere below we use the system of units with U = ks = !;• 



II. PHASE SEPARATION 



In order to take into account phase separation, we must admit such a possibility from the beginning, and then the 
corresponding phase probabilities must be defined in a self - consistent way. Assume that a phase separation has 
occurred so that our sample consists of regions occupied by two different phases numbered by the index = 1, 2 ; the 
phase regions being randomly distributed in space. Assign the index v = I to the superconducting phase, and ;/ = 2 , 
to the normal phase. For each phase distributicm. or phase configuration, the sample is nonuniform, which greatly 
complicates the consideration. However, assuming that each phase distribution is random, we may average observable 
quantities over these phase configurations. As a result, we come to a system described by a rcnormalizcd Hamiltonian 
that allows usual techniques supplemented with additional equations for phase probabilities and conditions of stability. 
The procedure of averaging over configurations has been expounded, in all necessary detail, in earlier papers [28-31] 
and discussed in a recent review [18]. Therefore, we shall not repeat it here, especially since the mathematical 
foundation is quite long and tedious even though the final result is rather simple. 

The averaging over phase configurations, in the case of a random mixture of two phases, leads to the definition of 
a renormalized Hamiltonian 



H=Hi®H2, (1) 

which is a direct sum of terms 

= w^h!;"^ + wIh:^\ (2) 

in which i?^*" is a kinetic - energy operator, including external fields, if any, iJ*"' is an interaction - energy 
operator. The renormalizing factor w^, is the geometric probability of the v - phase, that is, the ratio of the average 
volume occupied by the v - phase to the total volume of the sample. The Hamiltonian (2) is defined on the space 

y = yi(^y2, (3) 

which is a tensor product of Hilbert spaces, where y,^ is the space of states typical of the u - phase. When v = 1 
corresponds to a superconducting phase, and = 2 , to a normal phase, the state - spaces are chosen so that the so 
- called anomalous averages calculated with the states from 3^1 are not identically zero, while those calculated with 
the states from 3^2 are zero. In other words, the gaps related to the anom alous averages satisfy the condition 

Ai ^0, A2 = 0. (4) 
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This condition can be interpreted in several ways. For example, we can always posit that each Hilbert space of 
microstates, 3^^ , is restricted to functions which are invariant under a symmetry group yielding desired properties 
for the order parameters [32]. Then, the order parameters are just the gaps in (4), and the corresponding microstates 
can, for example, be chosen as the BCS wave functions [1,33,34] for the superconducting phase and the usual Slater 
determinants for the normal phase. Another way to interpret the choice of the order parameters (4) is to connect 
them, and the values of the anomalous averages, with the off - diagonal long - range order of reduced density matrices 
[35]. The largest eigenvalues of the latter are also known [35-37] to be directly related to this off - diagonal order. A 
more detailed discussion of these relations can be found in reviews [38,39]. Such relations between microscopic and 
macroscopic characteristics can be described in a general way by introducing the notion of order indices [40,41]. If 
is an n- particle reduced density matrix for a system of N particles, then the order indices are defined as 

lim 



Af— ►oo log A'' 

It has been shown [41] that the index q;„ satisfies the restriction 

< a„ < n (0 < n < A"). 

Different thermodynamic phases correspond to different sets of order indices, so that a v - phase is characterized by 
a set {ttra'^'} • Necessary and sufficient conditions that the phase with v = 1 is superconducting and the phase with 
u = 2 is normal are given by 

(1) (2) / fK\ 

Thus, there are several well developed methods for defining thermodynamic phases. What we still need to define is 
how to find the geometric probabilities Wv of the corresponding phases when considering a mixture described by the 
Hamiltonian (1). 

According to the general procedure [18], the phase probabilities are given by the minimization of the thermodynamic 
potential 

/ = In Tre-^^ (/3T ^ 1) (6) 

under the condition 

u;i+W2 = l, 0<w^<l, (7) 
where T is temperature. The normalization condition (7) can be taken into account explicitly by using the notation 

W\=W W-2. = \ — w. (8) 



Then w satisfies the equations 



dw ' dw^ 
The first of the equations in (9) is equivalent to 



-0, -L>Q. (9) 



0. (10) 



dw 

Substituting the hamiltonian (1) and introducing the notation 

^ {Htn, ^ (irr*>, (11) 

we obtain the equation 

2U2 + K2-K1 

for the probability of the superconducting phase. In case of thermodynamic phases with equal average densities, the 
probabilities Wi, coincide with the phase concentrations defined as the ratios of the number of particles in each phase 
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to the total number of particles. The average density of the superconducting and of the normal phase are practically 
the same. Therefore we may speak of (12) as of the superconducting phase concentration. 

The inequality in (9) shows when the phase separation in a sample is thermodynamically advantageous, as compared 
to a pure superconducting system. Taking account of (6), we find 

(0>)>«(lf)>- 

In accordance with (7) and (12), we have 

- 2C/i <Ki-K2< 2U2. (14) 
Substituting the Hamiltonian H , defined in (1) and (2), into the left - hand side of (13), we get 

= 2(f/. + t/.). 

Using this and noticing that the right - hand side of (13) is always non - negative, we obtain the inequality 

Ui + U2> 0. (15) 

Thus, phase separation is possible only when there are sufficiently strong repulsive interactions in the system for (15) 
to be satisfied. In other words, inequality (15) is a necessary condition for phase separation. This general result is in 
agreement with Hubbard - model calculations for the copper - oxide superconductors [42]. According to these it is 
necessary to include a sufficiently strong nearest - neighbour Coulomb repulsion in order to produce phase separation. 

A physical system will remain in a mixed state, with phase separation, as long as this is thermodynamically 
favourable compared with pure phases. The boundary, in the space of thermodynamic parameters, between mixed 
and pure phases is the set of points at which a necessary condition for phase stability is violated, - one of the conditions 
(13) or (14), for example. Experiment supports the usual assumption that transformation of a pure phase into a mixed 
phase begins with the appearance of nuclei of a competing phase within the pure phase. This process can be called 
nucleation and the point in the phase diagram at which this occurs, a nucleation point. At such a point, some 
thermodynamic and dynamic characteristics may display a singularity. This may occur, for example, in the density 
- density response function which is determined by the Fourier transform of the second - order density matrix. The 
behaviour of the thermodynamic and dynamic characteristics in the vicinity of such a phase transition is a separate 
problem which does not fall within the purview of the present paper. Rather, it is our object to delineate features 
which we regard as basic to any reasonable model of superconductors which exhibit mesoscopic phase separation. 

III. SPECIFICATION OF HAMILTONIAN 

Let us now specify the Hamiltonian, defined in (1) and (2), keeping in mind that Hi corresponds to the supercon- 
ducting phase; and H2 , to the normal phase. This implies that the corresponding order parameters and order indices 
must satisfy conditions (4) and (5), respectively. These conditions can also be formulated for anomalous averages 

defined through the field operators, Vi^s ( )) where s denotes spin and reR^. 
The anomalous averages 

(^is(r)Vis(r'))^0 (16) 
for the first representation are not zero, while those for the second representation vanish, 

(V'2.(r)V2.(r'))=0. (17) 

As is obvious, the conditions (16) and (17) are directly related to the corresponding properties of symmetry of 
microscopic states [4,18] or to the values (5) of the order indices [40,41]. 
The kinetic part of Hamiltonian (2) can be written in the form 

JJk^n ^ J Y^^lir ) [k47,7') - ij6{7 - /)] V..(rV 7 d 7', (18) 
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in which K^{r ^ r ) is a kinetic kernel including external scalar fields, if any, and /U is the chemical potential. Take 
for the interaction part of (2) a general expression 

HT'^\j 5^VL(^i)V'L(^2)K(ri, r2, rs, r4)V..'(^3)x 

ss' 

X il)us(ji)d ^1 (1^2 d^3 d r^4, (19) 

where V{. . .) is an interacting vertex including all effective interactions, repulsive and attractive, direct and indirect. 
For this moment, we do not need to specify these interactions. However, as an example, we can think of direct 
interaction as a repulsive Coulomb force taking screening into account, and that indirect interactions are mediated 
by an exchange of phonons or other bosons. 

In order not to complicate our discussion with the consequences of a variety of other known factors on the properties 
of superconductors, we limit ourselves to isotropic matter. Thus, we leave aside anisotropy effects and the related van 
Hove singularities [7]. For an isotropic system we can use the expansion 

V'.,(7) = -^^c.,(7)e^^^, (20) 
fe 

in which V is the total volume of the system. 

In what follows, to simplify the notation, we shall accept the convention 

cu(k) = Cs(k). (21) 

Then, we shall deal mainly with expressions corresponding to the superconducting phase, since an identical treatment 
can be given for the normal phase. The difference between the phases can be taken into account at the final stage by 
invoking conditions (16) and (17). 

Substituting (20) in (18), and using the property 



^ J n. hJ.h . 



we have 



k,k' ,s 



Ti{k, k ) - tJ-Skk' 



cl{k)cs{k ), (22) 



where the convention (21) is used and 

T^ik, k')=^j e-'^^K^{r, r ')e^^'^'rf r r ' (23) 
is the transport matrix. For (19) we get 

^^"* = ^ E ^^Sl,k2,k3,ki)cl{kl)cl,(k2)Cs'(k3)Cs(k4) (24) 



k\...kA 



with the vertex 



ri(fci, fc2, fcs, ki) = y"^i(r^i, r^2, r^3,^4)x 

X exp|-i fcir^i -i k2^2 +i ks^a +i ^4^4} d7 ^ d^ 2 dr 3 dl^ i . (25) 

In the following step to simplify the Hamiltonian, we use a convenient approximation which embodies a possible 
Cooper pairing. This could be done in several equivalent ways. We may construct an approximating Hamiltonian [43] 
acquiring the so - called Hartree - Fock - Bogolubov approximation. The latter is widely known and is successfully 
used in various applications, including such exotic ones as heated rotating nuclei [44-46]. We may follow a variational 
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approach as in the Blatt quasi - chemical equihbrium theory [47]. Or we may use the fundamental ansatz for the 
many - body wave function yielding what is called [38,39] the antisymmctrized geminal power. All these approaches 
are equivalent to one another [4,34,48,49]. To our mind, the main common idea lying behind these approaches can, 
in second - quantization language, be expressed by the operator approximation 



c\clc3C4 ~ C}C4(4C3) + {c[c4)c^C3 - {c[c4) {clca) - 

-c\c3{clc4) - (4c3)4c4 + {c\c3) {(^€4) + 
+ cj4(C3C4) + (c|4)C3C4 - {c\4;) {C3C4) , 



(26) 



where, for compactness, we write c„ instead of Cs„{k) ■ The right - hand side of (26) may be called the operator 
antisymmetrized geminal power. Averaging (26), we obtain the familiar approximation for the correlation function 



(c|4c3C4) = (4c4)(4c3) - (4c3)(4c4) + (4c^)(c3C4) 



of the Hartree - Fock - Bogolubov form. 

Applying (26) to (24), we present the latter as a sum 



in which the first term 



ki...k4 ss' 

X [4(fcl)c,(fc4)(4'(fc2)c.'(fc3)) -4(fcl)c«.(fc3)(cl,(fc2)c«(fc4)) 

has the normal Hartree - Fock form with 

Fjy {ki, /C2, k3, ki) = - ^v{ki, k2, ks, k4) + ^i^{k2, ki, k4, ks) 



The second term in (28), that is. 



^i"^ ~ 2Y ^ y^,ri(A;i, k2, ka, ki)"' 



ki...k4 



X 4(fcl)cUfc2)(c.'(A^3)c.(fc4)) +C,(fci)c,-(fc2)(cl,(fc3)4(fc4))] , 

contains anomalous averages and is, thus, responsible for superconductivity. The last term in (28) is 

-Bi = - ^2 y^r'i(A;i, k2, ka, fc4)x 



(27) 



(28) 



(29) 



(30) 



(31) 



{cl(kl)Cs'{k3)){cl,{k2)Cs{k4)) - {cl{kl)Cs{k4)){cl,{k2)Cs'{k3)) 
-{cl(kl)cl,(k2)){Cs'(k3)Cs(k4)) . 



(32) 



The Hamiltonian defined by (1), together with (2), (18) and (19), describes a system in which the total momentum 
and spin are conserved. Therefore 



{cl{k)Cs'{k )) = Skk'Sss'{cl(k)Cs{k)), 

{cl{k)cUt)) = S.kk'S-ss'{cl{k)cU- k))- 



(33) 
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If we denote the up and down spins as s =1, J, , then in (33) for s =t the notation —s means J, , and for s =J, , —s 
means t • For the normal average introduce the notation 



niik) = '^nu{k); nu{k) = {cl{k)cs{k)), 



(34) 



and for the anomalous average, 



ai(fc) = {c-s{- k)cs{k)). 



(35) 



The function nv{k) is a momentum distribution of particles. Since the transport matrix (23) and the vertex (25) 
do not depend on spin, we have 



n^{k) = 2n^s{k)- 

The normal term (29), with the use of (33) and (34), can be written as 

iff- = J2 Mi(k,k')cl(k)cs(k'), 

kk' ,s 



where Mi, is a mass operator with kernel 



M,(fc,fc') = ^E 



V 



(k, P,P,k)-^Ti. (k, P,k ,P) 



n.iP), 



By using (33) and (35), the superconducting term (31) becomes 



1 



jjsup^_J2r,{k,k,-P,P)x 
kk'p 

cl{k)cl,(k')ai(p) + cs{k)c-s(k') k (P) 



And the scalar term (32) simplifies to 



= ~^M^(k,k)nl(k)- 



-^Y.'^l{k,-k,-P,P) k {k)Ti(P)- 

kp 



(36) 



(37) 



(38) 



(39) 



(40) 



Notice that the necessity of an accurate and detailed analysis of all transformations presented here is dictated by 
our main goal which is to consider the role of all terms of the Hamiltonian not omitting any of them. As we have 
already shown and shall see in what follows, all these terms, repulsive as well as attractive, play a crucial role for 
superconductors with phase separation. 



IV. GAP EQUATION 

The sole approximation we have invoked so far now is the operator geminal power form (26), which is equivalent 
to the Hartree - Fock - Bogolubov approximation or to the antisymmetrized geminal power approach. The structure 
of the resulting Hamiltonian is still too complex to permit specific conclusions. Actually, a group theoretical analysis 
by Ozaki [50] of the Hartree - Fock - Bogolubov approximation, enumerated 73 possible ordered states, including 26 
superconducting states! The latter include non - Cooper pairing, when an electron pair has nonzero momentum, and 
various superconducting states coexisting with other nonsuperconducting orders. 

To make the problem tractable, we can resort to an approximation that restricts the space of microscopic states 
to a subspace satisfying some additional constraints [48,51,52]. In our case a natural such constraint is to consider 
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only those quantum states that conserve momenta, that is, when the momentum conservation occurs not only on the 
average, as in (33), but for all operator combinations: 

cl(k)cs(k') = Skk'cl(k)Cs(k), 

clik)cU{t) = S^kk'cl{k)cU- ~k). (41) 

Such a restricted space consists of BCS wave functions that are a particular kind of antisymmetric geminal power 
functions [34,53]. 

The restriction (41) makes it possible to greatly simplify all formulas. To this end, let us introduce the single - 
particle spectrum 

s4k) = n(k,k), (42) 

the diagonal mass - operator 

M,(k) = M,(k,k) (43) 

and an effective interaction 

J^(k,P) = T^(k,-k,-V,P). (44) 
Taking account of (41) and (42) - (44) we obtain for the kinetic term (22) 

iff- = ^ [ei{k) - /x] c\{k)cs(k), (45) 

for the normal term (37) 

= Y.Mi{k)c\{k)cs{k), (46) 

and for the superconducting term (39) 



Hr = ^12Mk,P) \cl{k)cU- k)a,{p)- 



kp,s 



+Cs(k)c-s{- k) ^1 (P) 



(47) 



The scalar term (40) becomes 

S. = -^^M,(fc)n,(fc)- :^^J,(fc,P)at(fc)a,(p). (48) 

k kp 

Finally, collecting all these terms, for the Hamiltonian of the superconducting phase, given in (2), we obtain the 
expression 

Hi=wiY,^Mk)cl(k)cs(k)+wlHr + wlBu (49) 

k s 

in which fl"™^ and Bi are defined by (47) and (48), and 

cj^{k) = e^(k)-fi + w^M^(k) (50) 

plays the role of an effective spectrum, renormalized, as compared to the single - particle spectrum (42), by the 
presence of the mass operator (43). 

The Hamiltonian (49) can be diagonalized by the Bogolubov canonical transformation 
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Cs{k) = u{k)as{k) +v{k) a-s (- k), 
in which both Cg and ttg satisfy the Fermi commutation relations. Diagonalization is achieved with 



(51) 



u{k) 



1 



1 + 



wi(fc) 



Eiik) 



v{k) 



1 



oJiik) 



Ei{k) 



leading to the Hamiltonian 



Hi=wi^Ei{k) k (k)as(k)+wiCi, 



where the quasiparticle spectrum 



contains the gap 



Et{k)^Aiik)+u;t{k) 



Ai(I) = -|i^Ji(fc,p)ai(P). 



For the scalar part of (52) we have 

Cu = J2 K(^) - ^''(^^ + ^AkWik) - ^M^(k)n^(k) . 

k 

With the diagonal hamiltonian (52) it is straightforward to calculate the momentum distribution (34), 



t7^ 1 t^l(fc) , , WiEi{k) 

ni{k) = 1 ^ tanh ■ 



Ei{k) 



2T 



and the anomalous averages (35), 



cri{k) = ^tanh- 



2Ei{k) 



2T 



So, for the gap (54) we obtain the equation 



^ ^ ,^ ^. Ai(p) , wiEi(p) 

Ai(fc) = -- r^2Z-^i(fe'^) ^ J tanh- ^ 



2Ei{p) 



2T 



(52) 

(53) 
(54) 

(55) 
(56) 
(57) 
(58) 



Let us emphasize that in the case of a superconductor with phase separation, with which we are dealing, all formulas 
obtained above involve the phase probabilities w^, in an intricate way. And it is very important where and how 
the latter enter into the expressions. Even though the technical approximations wc have adopted are standard, it 
is essential to accurately trace the role of the phase probabilities since it is these which distinguish the heterophase 
from pure superconductor. Also, we have to carefully take into account all interactions since, as has been shown, the 
presence of repulsive interactions is decisive for the phase separation itself. 

For the superconducting phase we need a nontrivial solution of (58). For the normal phase according to condition 
(17) instead of (56) we have 



n2{k) = 1 - tanh 



W20J2{k) 



2^ exp[w2iV2{k)/T] + 1 

and instead of (57) and (58) we have 

CT2(fc) = 0, A2(fc) = 0, 
with the effective spectrum, ui2{k) , given by (50). 



(59) 



(60) 
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V. INTERACTION POTENTIALS 



It follows from the previous Section that all characteristics of the system can be calculated provided we know the 
mass operator (43) and the effective interaction (44). These quantities are not independent of each other since they 
are both defined, by means of (38) and (44), through the vertex (25). If the real - space vertex entering into the 
hamiltonian (19), describes an interaction between particles, which is, as usual, invariant with respect to permutation 
of particles and space reflections, then the momentum - space vertex (25) has the properties [2]: 

ri^(fci, k2, ks, ki) = r^(fc2, ki, ki, fcs) = 

= r^(— ki,— k2,— ka,— ki) = ^u{ki,— ks,— k2, ki) = 

= ^uik3, ki, ki, k2)- 
With these properties, the symmetrized vertex (30) coincides with (25), 

Tv {ki, k2, ks, ki) =^iy{ki, k2, ks, ki)- 
It then follows from (38), (43) and (44) that 



M.(k) = ^^n,(p) 



V 
p 



Mk,k) - ^Mk,P) 



(61) 



Thus, the mass operator (43) is completely defined by the effective interaction (44) by the relation (61). 

As we have emphasized, we have to keep in the Hamiltonian (19) all interactions, direct and indirect. Therefore, 
the vertex (25) can be written as a sum = P^"' + F™*^ of direct and indirect interactions. Hence, the effective 
interaction (44) is a sum 

Mk, P) = Jf^ik, P) + JT'^ik, P) (62) 

of direct and indirect interactions. 

In principle, we could continue our analysis employing only the occurrence of the separation (62) of the effective 
interaction into two terms, without specifying their nature. For example, as interacting particles we coiild think of 
electrons or of atoms of '^He or of nucleons, taking for the real - space direct interactions the Coulomb, Lennard 
- Jones or Yukawa potentials, respectively. The indirect interaction might be assumed to be mediated by a boson - 
exchange mechanism involving phonons or excitons or something else. 

However, we prefer to be more concrete, so that in what follows, we shall think of particles as electrons whose 
direct interaction is described by a screened Coulomb potential and indirect interaction induced by phonon exchange. 
Generally, the screened Coulomb potential depends on the properties of a phase inside which electrons interact, since 
the screening is described by an inverse dielectric function reflecting the features of matter [54,55]. The simplest form 
of screening is obtained by using the Thomas - Fermi approximation which is equivalent to the static approximation for 
the Lindhard dielectric function [54,55]. Then, in real space the screened Coulomb interaction of particles immersed 
in the v - phase has the form 

Mr) = ^e-'=■'^ (63) 



in which e is the charge, k^^ is a screening radius with 



2 4 /3 y/' ( 1 

= — -n„ as 



where m is the electron mass and rii, , the electron density in the u - phase. In momentum space, (63) yields 
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Recall that the particular form (64) is not of great importance for our argument: for example, we could take for 
Vv{k) the more general expression (47r/fc^)£j7"'^(fc, wj;^^(fc)) , in which £y(fc,w) is a dielectric function and ujf,^'^{k) , a 
spectrum of elementary excitations in the v - phase. Also, for the dielectric function we could invoke a more refined 
approximation [55-60] . Here we have chosen the simple form (64) in order to illustrate how the screening can depend 
on the phase properties. In any case, whether Vi,{k) is given by (64) or by a more complicated expression, the part 
of (62) related to the direct interaction is 

Ju'^ik) = j $.(r)e-^(^-^)^d 7=v,(\k-P\). (65) 

For the effective indirect interaction we take one mediated by phonon exchange in the form [55] 

, P) = (66) 

where u^o is a characteristic phonon frequency, is the electron - phonon coupling taking account of screening, 



Q!l/ = — 



fc.f(l + K2/fc2^) 



Km) 



Z and M are an ion charge and mass, respectively, is the density of ions in the phase, and fcj/F ~ (STT^rii/)^/^ 
is the Fermi momentmn of an electron. 

The presentation (62) of the effective interaction as a sum of direct and indirect interactions leads to a similar 

decomposition 

M,(k) = Mf-(k) + Mt^'ik) (67) 
of the mass operator (61). The first term in (67) is related to (65) giving 

Mt'^{k) = n.v,{Q) - l^Y.n,(p)v,{\ k -P \), (68) 

p 

while the second term in (67) is expressed through (66) yielding 

M'J^'G) = -n}^ - 1 ^n.(p)4"'^(fc, p), (69) 

where 

p 

is the electron density in the u - phase. 

The chemical potential, as a function of temperature and density, is defined by the formula 

N = -{^)=Y,w,J{L{r)Mr))dr (71) 

for the total number of electrons in the system. Eq.(71), with notation (34) and (36), gives 

AT = ^ [^wini(p) + W2n2(p) ■ (72) 
p 

Using (70), we can reduce (72) to the equation 

N 

n = — = Willi + W2n2 (73) 



11 



connecting the average electron density in the system with the electron densities in each of the phases which the 
system is composed of and with the corresponding phase probabilities. Note that when the electron densities in both 
phases coincide with each other, that is when n,y ^ n , then (73) reduces to (7). 

Considering the gap equation (58) we pass to the thermodynamic limit with the standard replacement y — > 

/ (2;f).i • The main contribution to the resulting integral in the right - hand side of (58) comes from momenta close 
to the Fermi surface which is defined by the condition 



uJiikp) = 0. 

Keeping this in mind, we can rewrite the gap equation (58) as 

|2 



(74) 



Ai(P) 



— - vi (k) 
/^l(p) + uol(p) dp 



: tanh 



2T 



{2n) 



3 ' 



(75) 



where 



{k) = lim -3- /vi(| k-P \)dn(p) 
p—^kp 47r J 

"2 {k + kpf + AC? 



7re 
kkp 



In 



(fc - kpY + kI 



(76) 



is the screened Coulomb interaction averaged over spherical angles. 

A necessary condition that (75) has a nonzero solution for the gap Ai(p) is 



"1 



vi (fc) > 0, 



(77) 



which occurs in the vicinity of the Fermi surface so that uJi^k) < tof^ . Actually, inequality (77) has almost the same 
form as the BCS criterion for superconductivity [1]. However in our case the quantities entering into (77) depend on 
the superconducting phase probability. This dependence is essential and as we show below can dramatically change 
the characteristics of superconductors with phase separation making the existence of superconductivity possible even 
for values of the parameters which in a pure sample would imply the absence of superconductivity. 
Equation (75) can be simplified by considering the value of the gap at the Fermi surface, i.e. 



and by introducing the level density per spin. 



A= lim Ai(A;), 



dS{uj) 



(27r)3 



|V/cWl(fc) 



(78) 



where the integration is over the surface given by the equation wi ( fc ) = . Then (75) yields 



Jo 



X. 



■eff 



: tanh 



w 



where the effective coupling parameter is 

Xeff = wNi{0) 



2T 



dio 



■- vi (kp) 



(79) 



(80) 
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A^i(O) being the level density (78) at the Fermi surface defined by (74). 
The criterion (77) implies that for superconductivity the effective coupling parameter (80) is positive. For metals 

one has ki ~ kp and vi (kp) ~ ne'^/kp . Therefore in this case we have the estimate 



^1 (fcf ) 

l«ilV^? 



47r 



M' 



in which flp is the ion plasma frequency. The condition Ae// > means that wio < i^p ■ The latter inequality 
makes it possible to understand why structural instabilities and the related softening of the lattice, typical of high 
- temperature superconductors, favour the appearance of superconductivity. Lattice softening is associated with a 
decrease in the characteristic phonon frequency wio , this makes it easy to satisfy the condition u)iq < Qp and favours 
the onset of superconductivity. However, the latter occurs not in the whole volume of a sample, but only in those parts 
of it that are occupied by the superconducting phase. This is why lattice softening, enhancement of superconductivity 
and phase separation are phenomena that are intimately related to one another. This conclusion is supported as well 
by the consideration of s Hubbard - type model for which it has been found [42] that superconducting correlations 
are enhanced in the phase - separation regions. 



VI. CRITICAL TEMPERATURE 



The superconducting phase transition occurs at the critical temperature Tc where A = 0. . The equation for Tc 
follows from (79): 

f"'" 1 WOJ 

X,ff -tanh— du; = l. (81) 

Jo ^ ^^c 

Recall that here w{Tc) is the geometric probability of the superconducting phase, that is, the relative part of 
volume occupied by this phase. The upper limit in the integral (81) is the characteristic phonon frequency exhibiting 
softening because of structural fluctuations and the related phase separation [18,61]. The simultaneous occurrence 
of lattice softening and phase separation can be clearly observed by using Mossbauer spectroscopy [62,63]. In high - 
temperature superconductors these effects lead to anomalous saggings of the Mossbauer factor at critical temperature 
[12,13,64]. The softening of phonon frequencies at can also be observed by other experimental methods such as 
infrared, Raman and neutron scattering [11]. The softening of the characteristic frequency wio , as was shown in 
[61,63], can be expressed by the relation 

wio = w'^ojo, (82) 



in which loq is a characteristic phonon frequency of a pure superconductor, without phase separation, and the 
parameter ip measures the intensity of softening. Weak softening corresponds to (p = 1/2 , moderate, to (p = 1 , 
and the strong softening, to = 3/2 . 

For superconductors with phase separation the dependence of the critical temperature Tc on the superconducting 
- phase probability w is given by an intricate relation through (80) - (82). If one interprets the onset of supercon- 
ductivity as the condensation of Cooper pairs, then one can call w the concentration of superconducting condensate. 
No matter what it is called, the main point is that this quantity can be varied in high - temperature superconductors 
by varying their chemical structure, for example by doping. The experimentally measured dependence of Tc on w 
exhibits, for some high - temperature superconductors, a quite unusual nonmonotonic behaviour, as is reviewed in 
refs. [65,66]. For this reason it is especially interesting to study the dependence of Tc on iv . 

Formula (81) allows us to obtain two asymptotic expressions. One limit corresponds to the case 

Tc«—w^+'^, Ae//«1, (83) 

TT 

when 

Tc ~ 1.134«;^+*'a;oexp ( ---^ ) • (84) 

V ^eff J 

Note, that if we consider oJeff = w^'^'^ujo and Ae// as fitting parameters, then expression (84), as is known [67,68], 
describes the majority of low - temperature superconductors quite well. 
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Another limiting case, opposite to (83), is when 



TT 



Xeff » 1. 



Then (81) gives 



Te ~ -w'+'^XeffLJo. 



(85) 



(86) 



There is a temptation to treat (85) as the strong coupling limit. In doing this, we should be very cautions and not 
forget that Ae// , given by (80), is an cfFoctive coupling parameter. So, it may happen that Xeff is large, even 
if the bare coupling constant, A is not. If A » 1 , then by a standard argument Eliashberg equations imply that 
Tc ~ VX , although the dependence Tc ~ A , as is stated in [69], is also consistent with the strong coupling limit of 
these equations. 

The involvement of the superconducting phase concentration in the definition of the effective coupling parameter (80) 
makes the applicability of such simple asymptotic expressions, as (84) and (86), to high - temperature superconductors 
with phase separation quite limited. 

We therefore attempt an accurate analysis of the dependence of Tc on w . First, we have to remember that 
the level density (78) also depends on w through the effective spectrum (50) renormalized by the mass operator 
(67). Consider the isotropic case leaving aside the very interesting, but separate problem of van Hove singularities 
[7]. Then, differentiating the effective spectrum (50), we find 



lim VkUJiik) =e'p + wM'p, 



(87) 



where 



The level density (78) at the Fermi surface becomes 



^i(O) 



where we use the notation 



iV*(0) 



27r2| 



M'p 



|1 + 7w| ' 



7; 



(88) 



(89) 



To estimate the value of 

parabolic zone e'p ~ kp/m 



M'p 



take into account that n\ 

so that for the parameter 7 



k^p /Stt^ and ki 

I (89) we get 7 



kp , then M'p w 46^/377 . For a 
ae/ir'^aB , where is the average 
distance between electrons, and as is the Bohr radius. In good conductors ~ as , whence 7 ~ 0.1 . This means 
that in good metals the renormalization of the spectrum, due to the mass operator, is quite weak. 

In bad conductors with low electron density one has ^ Ub ■ Therefore 7 may equal 1 or more. This 
makes the role of the mass operator in renormalizing the electron spectrum and consequently in influencing the level 
density (88) very important. Such a situation is directly related to high - temperature superconductors which are, as 
is known, bad conductors having a low density of carriers. 

To distinguish the effects due to phase separation, we introduce the quantities 



A* = N^{0) 



2 ' 



(90) 



which do not depend on the phase probability w . The first quantity in (90) is the electron - phonon coupling constant 
for a pure system without phase separation. The second quantity in (90) is the average intensity of screened Coulomb 
interaction. 

With notation (90), for the effective coupling parameter (80) we obtain 



A, 



eff 



X* - IJL*W^^ 

vfif-^l + 7w| ■ 



(91) 
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The criterion of superconductivity means that Ae// > , which leads because of (91) to 



A* - ^Ji'w'^'f > 0. 

This inequaUty, in the case of pure superconductor, when w = 1 , reduces to the standard condition A* > ji* 
usually valid for low - temperature superconductors. When the phase separation occurs in a superconductor, then 
the condition Ag// > is easier to satisfy. Really, when < w < 1 , then the inequality A* > ji*w'^'^ can 
be true even if A* < /u* . In this way, phase separation favours the appearance of superconductivity. A sample 
with phase separation can become superconductive even if a similar sample, without phase separation, cannot have 
superconducting properties. 

The superconducting critical temperature is given by (81) which becomes 



r ° 1 wu 

/ -tanh— da;=l. (92) 



As we see, it is not easy to analyse precisely the influence of phase separation on the critical temperature, that is, to 
solve explicitly for the dependence of Tc on the phase probability w . To slightly simplify the situation, we may 

notice that since e'p « kF/m, Ni{0) « rnkF/^ir'^ and Vi [kp) ~ ne'^/kp , then /i* « 7 . Therefore, for simplicity, 

we put 7 = /i* . 

Of course, such a slight simplification docs not help much, and to proceed further in defining the dependence of 
Tc on w we have to resort to numerical analysis of (92). To this end, it is convenient to introduce the dimensionless 
critical temperature 



tc — • 

Wo 



Now we can reorganize (92) in the form 



— ^ / tanh — — X — = 1. (93) 

w^'P-i{l + ij,*w) Jo \ 2tc J X ^ ' 

Prom (93) we can immediately conclude that the behaviour of tc = tc{w) can be, in general, nonmonotonic, since 
there are two points at which tc{w) tends to zero. The flrst point corresponds to w — > . Then from (93) we obtain 

l,_,A-^_ 

2 l + n*w ^ ' 

The second case, when ^ , is when w — > (A*//^*)^/^'^ , then 

Recall that the dependence of tc on w is interesting to analyse because the superconducting phase concentration 
w can be measured and can be varied in experiments by changing the chemical structure of materials, for example, 
by doping [6,11,65,66]. 

We made a detailed analysis of the function tciw) by solving the equation (93) numerically. Graphs of the resulting 
functions are presented in figs. 1-12. Wc did not try to fit any particular experimental situation. Rather we wanted to 
understand the whole picture as to how the behaviour of tc(w) changes qualitatively with the change of parameters 
ip^ix* and A* . 

We were pleasantly surprised by the wide variety of curves which were obtained. It is apparent from the accompa- 
nying graphs that by the choice of the corresponding parameters it will be possible to obtain a reasonably good fit to 
any experimental curve. 

Figs. 1-4 correspond to the case of weak softening; figs. 5-8, to that of moderate softening; and figs. 9-12, to the case 
of strong softening. With the increase of the Coulomb parameter the function tc(w) really becomes nonmonotonic. 

The behaviour of tc(w) plotted in figs. 3, 4, 7 and 8 has a striking similarity to the corresponding experimental curves 
for high - temperature cupratc superconductors (sec [6,11,65,70] and references therein). The coincidence becomes 

practically complete if we redraw our figures in the relative coordinates T= tc/t^""^, a = w/w"^""^ , as in [65,66], 
where the point (i™«»:^ denotes the point of a maximum of the considered curve. 
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It is clear from the proceeding analysis, as illustrated in our graphs, that, for fixed (p, /i and A , the value of the 
reduced critical temperature, Tc , depends crucially on the parameter w . The operational meaning of this parameter 
is rather different according as the the system is 1) stable, or 2) metastable. 

Notice that while w occurs explicitly on the left hand side of equation (12) it also occurs implicitly in a complicated 
manner in the right hand side as a result of the renormalization of the Hamiltonian (1). Thus for a stable system 
satisfying conditions (9) and (13), we can think of w as determined self - consistently by equation (12) once all 
characteristics of the system, - such as chemical composition, particle masses, interaction potentials, temperature, 
density and pressure - are given. Note that all of these characteristics are necessary for w to be uniquely determined. 
In particular, the total role of interactions must be taken into account. In contrast to this, following common practice, 
in our model the gap equation (79) and consequently equation (93) for the critical temperature, take into account 
only interactions specified on the Fermi surface. Essential to the argument of the present paper is the conviction 
that it would not be reasonable to try to treat superconductivity in heterophase system with a model in which w is 
determined by parameters defined merely on the Fermi surface. For instance, parameters characterizing the ground 
state would be indispensable. The equations for phase probabilities always contain more characteristic constants than 
the equations for an order parameter [18]. It is therefore permissible to contemplate the possibility of holding some 
parameters, such as A* and /x* , fixed while w varies as a result of changing other parameters such as the chemical 
composition. 

On the other hand, for a metastable system for which (13) does not hold even though thermal and mechanical 
stability are preserved, the fraction w is not necessarily determined by equation (12), but might be arbitrary 
depending on the preparation of the sample. This possibility should not be overlooked since many high - 
superconductors are metastable. 

VII. CONCLUSION 

We have developed an approach, for describing superconductors, taking into account three mutually interrelated 
factors: (i) repulsive interactions, (ii) lattice softening and (iii) phase separation. We have deliberately limited 
ourselves to the use of only commonly accepted approximations. This was in order to emphasize that our results are 
not artifacts of some technical tricks, but the direct consequences of the physical reasons offered. The main results 
can be summarized as follows: 

(i) A necessary condition for phase separation in a superconductor is the presence of repulsive interactions. 

(ii) Phase separation favours superconductivity making it possible in a heterophase sample even if it were impossible 
in a pure sample. 

(iii) The superconducting critical temperature as a function of the relative concentration of the superconducting 
phase can display the nonmonotonic behaviour typical of high - temperature cuprate superconductors. 

It should be a straightforward matter to adapt the basic approach of this paper to other models of superconductors 
by using alternative approximation methods. 
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Figure Captions 



Fig.l. 

The superconducting critical temperature as a function of the superconducting phase concentration for the param- 
eters ifi = 0.5, fj,* =0.1 and A* = 1 (lower curve), A* = 5 (middle curve) and A* = 10 (upper curve). 



Fig.2. 

The same as in fig.l, but for the parameters (fi = 0.5, fj,* = 1 , and A* = 1 (lower curve), A* = 5 (middle curve) 
and A* = 10 (upper curve). 



Fig.3. 

The same as in fig.l for the parameters if = 0.5, /x* = 5 , and A* = 1 (hardly visible lower curve). A* = 5 
(middle curve) and A* = 10 (upper curve). 



Fig.4. 

The same as in fig.l for the parameters (p = 0.5, /x* = 10 , and A* = 5 (lower curve) and A* = 10 (upper curve). 
The curve corresponding to A* = 1 in this case is invisible. 
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Fig.5. 

The same as in fig.l for the parameters ip = 1, /i* = 0.1 , and A* = 1 (lower curve), A* = 5 (middle curve) and 
A* = 10 (upper curve). 



Fig.6. 

The same as in fig.l for the parameters if = 1, yu* = 1 , and A* = 1 (lower curve), A* = 5 (middle curve) and 
A* = 10 (upper curve). 



Fig.7. 

The same for the parameters f = 1, /x* = 5 , and A* = 1 (lower curve). A* = 5 (middle curve) and A* = 10 
(upper curve). 



Fig.8. 

The same for the parameters <fi = 1, /U* = 10 , and A* = 1 (lower curve). A* = 5 (middle curve) and A* = 10 
(upper curve). 



Fig.9. 

The same for the parameters (p = 1.5, /j.* = 0.1 , and A* = 1 (lower curve). A* = 5 (middle curve) and A* — 10 
(upper curve). 



Fig.lO. 

The same for the parameters ip = 1.5, /U* = 1 , and A* = 1 (lower curve). A* = 5 (middle curve) and A* = 10 
(upper curve). 



Fig. 11. 

The same for the parameters (p = 1.5, /x* = 5 , and A* = 1 (lower curve). A* = 5 (middle curve) and A* = 10 
(upper curve). 



Fig.l2. 

The same for the parameters ip = 1.5, /U* = 10 , and A* = 1 (lower curve). A* = 5 (middle curve) and A* = 10 
(upper curve). 
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